Numerical studies on chaoticity of a classical hard-wall billiard with openings 
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Using fractal analysis, we investigate how the size of openings affects the chaotic behavior of a 
classical closed billiard when two openings are made on the boundary of the billiard. This kind of 
open billiards retains chaotic properties of original closed billiards when openings are small compared 
to the size of the billiard. We calculate the fractal dimension using an one-dimensional subset of all 
possible initial conditions that will produce trajectories of a particle injected from an opening for 
the billiard, and then observe how the opening size can change the classical chaotic properties of 
this open billiard. 
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Two-dimensional (2D) hard-wall billiard systems have 
been a popular subject for studying the dynamics of 
chaotic systems 0. Theoretically, 2D billiard systems 
can be used to study manifestations of classical chaos in 
semiclassical and quantum mechanics. The classical dy- 
namics of a billiard system shows three distinct types of 
behavior: the system is either integrable (regular behav- 
ior) or non-integrable (either soft chaos, characterized by 
mixed phase spaces that have both regular and chaotic 
regions, or hard chaos, characterized by ergodicity and 
mixing) M. In experiments, on the other hand, the bil- 
liard systems can be used as models to explain fluctuat- 
ing behavior of magnetoconductance through 2D semi- 
conductor heterostructures in the ballistic regime (H f|. 
In such an experiment, the structure have the shape of 
a 2D billiard, of which the dynamic properties of the 
billiard are already known (for example, the circle as a 
regular system, or the stadium as a system with hard 
chaos). To measure the conductance through the struc- 
ture, leads (from here on, we will call them openings) 
are attached to the structure; hence the whole structure 
becomes an open billiard. One usually studies transmis- 
sion properties through these open billiards, and links re- 
sults to the classical chaotic properties of original closed 
billiards, because the properties of the closed billiard 
are still observed indirectly when the size of the open- 
ings is small compared to the size of the billiard (see 
Refs. H ^, 0). Therefore, to better relate the open bil- 
liard to the original closed billiard, quantitative analysis 
on effects of openings in classical dynamics is necessary. 

In this paper, we numerically calculate a quantity that 
represents the chaoticity of classical open billiards, and 
see how the size of the openings affects the overall clas- 
sical dynamics. In Ref. the fractal dimension was 
obtained by looking at a two-dimensional set of initial 
conditions that will produce trajectories injected from an 
opening of the Sinai billiard with two openings. [A set of 
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initial conditions in four-dimensional phase space can be 
reduced to a two-dimensional set, because (1) energy of 
the particle does not affect the trajectory of a particle, 
and (2) initial locations of the particle can be restricted 
to an one-dimensional (ID) curve.] Here, however, we 
will use an ID subset of all possible initial conditions for 
calculations of the fractal dimension to make numerical 
calculations simple. The incident angle of the particle 
injected from the center of one opening represents an 
ID subset of all possible initial conditions, which give 
rise to some possible trajectories of injected particles. In 
Ref. ||, graphs of the "exit opening" (an opening from 
which the particle exits) vs the incident angle from the 
center of an opening were used. There are "geometri- 
cal channels" [jlOj [llj inside a range of incident angles. 
When there are two openings, we can define transmis- 
sion windows that represent initial angles for transmit- 
ting particles. The boundaries of these windows will form 
a fractal^], [|. In this paper, instead of calculating the 
fractal dimensions of boundaries of these transmission 
windows, we use the graph of the number of collisions vs 
the incident angle, and calculate the fractal dimension of 
the set of singularities, at which the number of collisions 
changes abruptly, inside the set of all possible incident 
angles (see Ref. fl^l). The fractal dimensions from these 
two graphs will be eventually the same. 

The billiard used here for numerical calculations is the 
circular billiard with a straight cut (the cut- circle bil- 
liard) with two openings (see Fig. Q). There are five 
parameters: (1) the width W measured in the direction 
perpendicular to the cut, (2) the radius R, (3) the angular 
width A of the openings, which represents the opening 
size, (4) the orientation angle of the cut relative to the 
first opening, and (5) the position of the second opening 
relative to the first opening as measured by the angle 7. 
We scale the width W by R so w = W/R, and thereby 
reduce the number of independent parameters to four: 
w, A, fl, and 7. For all subsequent discussions, we set 
O = 135° and 7 = 270°. For these values, w cannot be 
less than [1 — cos(45°)] ~ 0.293, and A has an upper limit 




when A = 30 (deg) 



FIG. 1: The geometry of the cut-circle billiard. The billiard 
has two openings, I and II. The size of the cut is given by the 
width W (we define w as W/R to represent the relative size 
of the cut). The size of two openings is A, and the position 
of the cut and the opening II, Q, and 7, are measured from 
the position of the opening I. A particle is injected from the 
center of the opening I with an angle cj>. 
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FIG. 2: Graphs of the number of collisions vs the incident 
angle <f). Here the number of collisions is scaled for each w 
to fit in graphs with the same sizes, because we are only in- 
terested in positions of singularities. Five cases with different 
cut sizes are shown : w — 0.5, w — 0.75, w = 1.04, w = 1.5, 
and w — 1.71. 
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2 [cos _1 (l -w)- tt/4] when 0.293 < w < 1 
7r/2 when 1 < w < 2 

(1) 

in radian. When the billiard is closed (A = 0), it has 
been proven that the phase space is mixed (soft chaos) 
when < w < 1, and that the phase space is fully chaotic 
(hard chaos) when 1 < w < 2[^3). The system is inte- 
grable when w = 1 and 2. This billiard has a property 
of showing all three types of dynamic behaviors just by 
changing a parameter w, and this is the reason why the 
cut-circle billiard is a good example for analysis like this. 

The particle is injected with an incident angle <j> 
(—it/2 < (f> < 7r/2). One can easily calculate several 
quantities by numerically following these trajectories, but 
we focus only on the number of collisions with the wall 
before the exit with respect to <fi. In other words, we find 
a mapping from the incident angle in opening I to the 
number of collisions before the exit. In Fig. ||, we show 
how the number of collisions changes with <p for five dif- 
ferent w values (w = 0.5, w = 0.75, w = 1.04, w = 1.5, 
and w = 1.71) when A = 30°. In these calculations, there 
are finite number of singularities when w = 0.5, but, in 
other cases, there are regions where singular points are 
closely packed together showing infinitely fine structures. 
We observe that there tends to be more singularities as 
w increases. That is because (1) when w is bigger than 
one, the original closed billiard shows hard chaos, which 
means global chaos, and (2) the opening size is fixed here 
for all u>'s even though the overall size of the billiard in- 
creases as w increases, as a result, making the relative 



size of the openings smaller. 

We calculate the fractal dimension df for sets of singu- 
larities for various cases, using a simplified box-counting 
algorithm. For Np uniformly distributed incident angles, 
we numerically find the number of collisions for each in- 
cident angle representing an initial condition. For sev- 
eral Np values, we find the number of singularities N$ 
(assuming that there only exists a singularity when the 
number of collisions changes for two successive incident 
angles). Then the fractal dimension df is defined by 
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logip N s 
log 10 Np ' 



(2) 



which is the slope of the graph of log 10 Ns vs log 10 Np. 
We use the ordinary least-square fit to find the slope us- 
ing points in the graph with limited Np. Because the 
numerical accuracy for numbers is limited, Np have an 
upper limit. Moreover the numerically given error of 
an initial condition will amplify as the particle travels 
through the billiard and hits the boundary while travel- 
ing. Here we use the double-precision numbers with 15 
significant digits, and then Np should be less than 10 15 , 
especially when the number of collisions is big. In most 
cases in our example, Np up to 10 7 were used to obtain 
d/'s with two significant digits. In some cases (for ex- 
ample, when w — 1.71 and A > 59°), however, iVp's up 
to 10 s were used to find d/'s, because the slopes didn't 
converge sufficiently when Np is less than 10 7 . 

In Fig. ||, we fix the width w at five different values used 
in Fig. and vary the opening size A up to A max for 
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FIG. 3: Graphs of the fractal dimension df vs the opening 
size A. For each w (w = 0.5, w — 0.75, w — 1.04, w = 1.5, 
and w = 1.71), A is varied up to a maximum value A max . We 
can observe that df does not always decreases monotonically 
as the opening size A increases. 



each w. An interesting feature of these curves is that they 
are not monotonically decreasing as A increases. Small 
peaks appear even after df reaches zero at a certain A 
value. For example, the c?/-curve for w — 1.04 reaches 
zero near A ~ 35°, but there is a small peak near A ~ 
79°. The reason for this phenomenon is that there are 
more possible trajectories when the opening is bigger. In 
our method for the circular billiard, the launching point 
of the particle gets closer to the center of the billiard as 
the opening size increases. 

In Fig. ^, two trajectories that doesn't exist when 
openings are smaller are shown. Since these trajecto- 
ries have the higher numbers of collisions and there is 
sensitive dependence on initial conditions near them, df 
becomes non-zero. In Fig. (§a), the billiard (w = 1.04, 
A = 79°) lost most of the original shape of the closed bil- 
liard, but trajectories near the one shown (4> — 84.44°) 
still make up an infinitely fine structure. In Fig. ||(b), 
the billiard (w — 1.71, A = 59°) has a trajectory 
(<f> = 88.40°) that gets close to the 3-bounce closed orbit 
in the circular billiard. (The closed orbits in the circular 
billiard are stable, but only marg inally(l4| |l§.) Even 
though only few bounces are shown, this trajectory ac- 
tually bounces more than two thousand times as the tri- 
angular shape rotates clockwise little by little for each 
rotation, before the particle finally exits through cither 
one opening. A high peak near <f> — 59° for w = 1.71 in 
Fig. H is due to trajectories near this one. 

In Fig. ||, we calculate d/'s as the width w varies from 
0.5 to 2 with the step size 0.01, for several opening sizes 
(A = 0.5°, A = 1°, A = 5°, A = 10°, A = 20°, A = 30°, 
A = 40°, and A = 50°). We can compare graphs in two 
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FIG. 4: Two trajectories showing why the function of df 
with respect to A does not always decreases monotonically. 
As the opening size A increases, new initial conditions for 
trajectories, which can be very sensitive to the incident angle 
4>, may appear. For example, (a) w — 1.04, A = 79°, and 
$ = 84.44°. (b) w = 1.71, A = 59°, and <j> = 88.40°. 



different regions (found in the closed cut-circle billiard): 
< w < 1 (soft chaos) and 1 < w < 2 (hard chaos). 
When «; = 1 and w — 2 (integrable cases), df is zero 
in our calculation. As the opening size A approaches 
zero (see graphs for A = 0.5° and A = 1°), we observe 
that the fractal dimension df approaches one in both 
regions in our calculation. When the opening size A is 
not big (see graphs for A = 5°, A = 10°, and A = 
20°), we observe that behavior in two regions is clearly 
distinct. In the region < w < 1, there are fluctuations, 
which comes from the mixed phase space structures of the 
billiard, and in the region 1 < w < 2, graphs are smooth 
because the phase spaces of the billiard have no structure 
due to ergodicity. On the other hand, when the opening 
size gets bigger (see graphs for A = 30°, A = 40°, and 
A = 50°), the distinction between two regions, observed 
in cases with smaller openings, starts to disappear; there 
are fluctuations in both regions. Ergodicity in closed cut- 
circle billiards with hard chaos is no longer an important 
factor in the dynamics of open billiards when openings 
are big, as expected. 
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FIG. 5: Graphs of the fractal dimension df vs w for 8 dif- 
ferent opening sizes, A = 0.5° (o), A = 1° (□), A = 5° (O), 
A = 10° (A), A = 20° (v), A = 30° (+), A = 40° (o), 
and A = 50° (□). When A = 5° and A = 10°, the behavior 
of the graph is clearly distinct for two regions: hard chaos 
(1 < w < 2) and soft chaos (0 < w < 1). But when A > 20°, 
the distinction between two regions disappears. 



We have so far observed numerical results for the cut- 
circle billiard with two openings. A method that cal- 
culates the fractal dimension using an ID subset of all 
possible initial conditions has been used to determine 
effects of openings on chaotic closed billiard quantita- 
tively. From numerical results, we have found that the 
fractal dimension as a function of the opening size is not 
a monotonically decreasing function, because there are 



more trajectories, which can be sensitive to initial con- 
ditions, for bigger openings. These new trajectories can 
make up infinitely fine structures, and subsequently can 
cause the fractal dimensions to increase with A. We also 
have found that, as the opening size gets bigger, the dis- 
tinction caused by soft chaos and hard chaos of the orig- 
inal closed billiard starts to disappear, because the open 
billiard that are constructed by introducing openings to 
the closed billiard starts to lose the original shape of the 
closed billiard as the size of openings gets bigger. 

There are several points worth mentioning. First, if 
we use the 2D full set of initial conditions to calculate 
the fractal dimensions as was done in Ref. df will 
be in the range of 1 < df < 2, but the overall behavior 
of df is expected to be similar to that of the results ob- 
tained here. Second, we can use a mapping from initial 
conditions to other quantities like dwell times (or travel 
distances) instead of the exit openings and the number 
of collisions. Even though the dwell time will be more 
expensive to calculate numerically, the set of singulari- 
ties from the dwell time will be almost the same as those 
from the number of collisions. Finally, it has been found 
to be important and necessary to understand relations 
between classical dynamics and quantum (or scmiclassi- 
cal) dynamics of billiards |l6 , 17, 1§[, because the ballistic 
scattering is usually realized in microscopic scales. Then 
it will be interesting to ask how this kind of fractal di- 
mensions can be related to the quantum and semiclassical 
dynamics of the same kinds of open billiards. 
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